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' A novel approach for studying phase transitions in systems with quantum degrees of freedom 

, is discussed. Starting from the microscopic hamiltonian of a quantum model, we first derive a 

' set of exact differential equations for the free energy and the correlation functions describing the 

^SJ ' effects of fluctuations on the thermodynamics of the system. These equations reproduce the full 

r I , renormalization group structure in the neighborhood of a critical point keeping, at the same time, full 

^ l' information on the non universal properties of the model. As a concrete application we investigate 

I the phase diagram of a Heisenberg antiferromagnet in a staggered external magnetic field. At long 

wavelengths the known relationship to the Quantum Non Linear a Model naturally emerges from 
our approach. By representing the two point function in an approximate analytical form, we obtain a 
closed partial differential equation which is then solved numerically. The results in three dimensions 
are in good agreement with available Quantum Monte Carlo simulations and series expansions. 
More refined approximations to the general framework presented here and few applications to other 
models are briefly discussed. 
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I . I. INTRODUCTION 

Q . In recent years several works have appeared in the literature on the quantum Heisenberg antiferromagnet, with 

Q ' particular attention to the two dimensional case, because of its deep connections to the physics of high temperature 
superconductors. Although the corresponding classical model was studied by means of several technkjues ranging 

^ ' from high temperature expansionstl to Monte Carlo (MC) simulationsH and Renormalization Group (RGp ap^ppaches, 

^ \ accurate numerical results on the phase diagram of the quantum model have become available only recentlyH'u. Since 

■ the pioneering works basedjOn spin- wave theoryO and on high-temperature exnaaisionsB, a better understanding of the 
' behavior in two dimensionsQ and of the critical properties in three dimensionsffl has been provided by more and more 

■ efficient MC methods. Nevertheless, several problems are still open, especially in the low and intermediate temperature 
\l ' region. This regime has been successfully investigated by means of the RG, after an appropriate mapping onto the 

, nonlinear sigma model [NLaMp. The satisfactory agreement between RG and experiraeptal results, however, seems 
to hold only for spin 1/2 and quantitative deviations between theory and experimentsiij appear for larger values of 
S. Both Monte Carlo simulationaO and series expansionsEJ suggest that the breakdown of the A^LcrAf description of a 
Quantum Antiferromagnet is due to the intrinsic limitations of the methods based on effective actions. In fact, these 
\ approaches are devised for describing the long wavelength behavior of the model and can be justified in a regime 
I 1 characterized by strong correlations, while the experimental investigations span a wide temperature range, including 
'"O ' regions where local antiferromagnetic order is weak. The same conclusions are also found by use of a new semiclassical 
Ph , approach in which the problem of the quantum Heisenhesg antiferromagnet is reduced to a classical one at an effective 
O ■ temperature taking into account quantum fluctuationstHl. 

Besides the problems concerning the Heisenberg antiferromagnets in bipartite lattices, several other questions are 
still open in the framework of phase transitions in quantum systems. Generally, accurate analytical approaches at 
finite temperature are not available and quantum simulations are often hindered by the "sign problem" . For instance, 
frustrated antiferromagnets or fermionic systems, like the celebrated Hubbard model, are still beyond the capabilities 
of today's simulation algorithms. Interestingly, this class of models might be investigated by a straightforward 
generalization of the method we are going to present here, at least in order to estimate the magnetic phase boundaries 
or the occurrence of phase separation. 

In this paper, rWje present a new approach to the phase transitions in quantum systems. A preliminary account can 
be found in RefJi3. We develop a general formalism which can be used in a wide range of temperatures both above 
and below a possible critical point. The method is then applied to the specific case of Heisenberg antiferromagnets, 
for different values of the spin and of the spatial dimensions. This approach, which we called QHRT (Quantum 
Hierarchical Reference Theory) is a quantum extension of the classical HRT which was developed in the context of 
liquid-vapor transitions and successively applied to several system, like Ising models and binary mixturealil. Our 
method relies on the fact that, for a large class of quantum systems, an exact mapping exists between the quantum 
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partition function and a fictitious classical one to which HRT can be applied. The main advantage of HRT, also 
shared by QHRT, is its microscopic character, which allows to obtain both universal and non universal quantities of 
the system under consideration. Non classical critical exponents together with a precise determination of the critical 
temperature and a correct description of the first order phase transitions are in fact obtained within the QHRT 
formalism. As already mentioned, the theory is very general and can in principle be applied to several physical 
systems, bosons, fermions and spins. 

Following the idea of the RG approach, an evolution equation is derived, starting from the perturbative expansion of 
the Helmholtz free energy. Fluctuations over progressively larger lengthscales are taken into account by varying a cut- 
off parameter thereby connecting the system with no fluctuations (mean-field approximation) to the fully interacting 
one. The equation obtained in this way is formally exact, although it is not written in a closed form. In order to 
get a closed equation, the knowledge of the structure of the dynamic correlation functions of the model as a function 
of frequency and wavevector is needed. This represents the only approximation present in the QHRT approach and 
requires a thorough analysis of its physical implications. Near a critical point our partial differential equation acquires 
a universal form and reproduces the RG structure for a three component order parameter with critical exponents 
exact to first order in e = 4 — d. Also as T ^ 0, the QHRT approach is able to predict the known behavior originally 
derived by a RG analysis of the NLaM. We solved the equation for several values of the spin, both in two and three 
spatial dimensions. In d = 3 we obtained an excellent matching with the available numerical simulation data, also 
from a quantitative point of view. 

The paper is organized as follows: In Section || we give the derivation of the theory. In Section IH the theory is 
applied to the Heisenberg antiferromagnets, the classical Heisenberg model being obtained as the limiting S ^ oo 
case. In Section |^ an explicit closure is given and the implications of this choice are discussed. In Section ^ an 
analysis of QHRT equation in proximity of a phase transition is given. In Section [Vl| th e explicit results, obtained 
via numerical integration of the evolution equation, are presented. Finally in Section |VH| conclusions are drawn and 
further possible applications of our theory are proposed. 



II. DERIVATION OF QHRT EQUATIONS 

In this Section we develop a general formalism for the description of critical phenomena in quantum systems. Our goal 
is to implement the basic ideas of renormalization group to the microscopic, many body hamiltonian of an interacting 
quantum model. The class of hamiltonians we are going to consider, in arbitrary dimension d, has the structure: 

H = Hr + V = Hr + ^ J dxdyp{x)w{x - y)p{y) (1) 

where Hr is a reference part whose properties are supposed known and the coordinates x and y span the d-dimensional 
space. V represents an interaction term which couples the (local) obscrvables p(r) at two different points via the 
translationally invariant two body potential w{x — y). In the following we will assume that: i) the local operators 
p{r) commute at different points and ii) the interaction w{r) can be Fourier transformed. Instead, no limitation is 
imposed on the form of Hr, which is not restricted to describe free particles. This class of hamiltonians is indeed 
rather general and includes several systems of widespread interest, like quantum magnets (where p{r) represents a 
spin variable) but also models of interacting fermions, both on and off lattice (where p{r) is the local density operator) 
and even systems with electron-phonon coupling, like the Hubbard-Holstein model. Generalizations to other cases, 
for instance models with many particle interactions, are possible but will not be considered here. 

As a first step, we will write the partition function of the quantum model in a form identical to that of a classical 
system. The latter will be then studied by means of the powerful methods already available in the framework of 
classical critical phenomena. This quantum to classical mapping is most easily carried out by expanding the grand 
canonical partition function S as a power series of w{r\— The method we use closely follows the standard treatment 
of the path integral formulation of quantum mechanicalj, although we carefully avoid to introduce any smoothness 
assumption on the classical fields. Therefore, our treatment is not limited to long wavelength or low energy and the 
equations we are going to obtain are formally exact in the whole phase diagram of the model. 

The grand partition function of the quantum system, S = Tr exjp{—(3H), is first written in terms of the (imaginary 
time) evolution operatorEJ as: 

Tr exp(-/3i7) - Tr {exp{-pHR)U{(3)} (2) 
where U (t) is the solution of the differential equation: 
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(3) 



with U{Q) — 1 as initial condition. V{t) is the time dependent operator: 

V{t) = eMHRt] V eM~HRt] (4) 
Equation (||) can be formally solved by iteration with the result: 

U{t)^y^^—^ dh... dt„Tt[Vih)...Vitn)] (5) 
„=o -^0 -^0 

where Tt is the usual time ordering operator. By inserting in (^) the explicit form of V{t) 

V{t) = i y dxdy J dt'p{x, t)w{x - y)S{t - t')p{y, t') (6) 

equation can be written as a power series of the {d + l)-dimensional interaction 

X') = w{x - x')S{t ~ t') (7) 

Here, the variable X identifies the pair {x, t) where x spans the d-dimensional coordinate space and t belongs to the 
interval (0, /?). The explicit expression for the grand partition function S is: 

= S«EVt / dXidX[...dX,,dX'^pR{Xi,...X'J<f>{Xi,X[)...<j>{Xr,,X:^) (8) 

n=0 ^' •' 

Here and in the following R labels the reference system defined by Hr. The imaginary time dynamic correlation 
functions appearing in Eq. pi^i, ••■^n) are defined, for a general hamiltonian H, by: 

...a;„,t„) = ti).../5(a;„, i„)) = ^^Tr je"'^-^ Tt[p(xi, ti).../5(a;„, t„)]| (9) 

where p{x, t) is the local density operator in imaginary time: 

p{x, t) = exp[Ht] p{x) exp[-Ht\ (10) 



The correlation functions in the perturbative expansion are given by Eqs. (|9|,|10|) specialized to the reference 
hamiltonian Hr. Since p{r) commutes on different sites, the functions are symmetric under the permutation of 
the labels and can be viewed as the correlation functions of a hypothetical (d + 1) dimensional classical reference 
system. The quantum nature of the problem is then contained in the structure of the reference system through "E.^ 
and pii{Xi, ...Xn). It is interesting to notice that the n-particle correlation functions p{Xi, ...X„) associated to a given 
quantum hamiltonian H can be equivalently obtained as successive functional derivatives of a generalized partition 
function S[J] with respect to a space and time dependent external field J{x, t) coupled to the quantum operator p{x). 
More precisely, the generating functional S[J] is still defined as the trace of an imaginary time evolution operator: 

S[J] =Tr{cxp[-/37?];7(/3)} (11) 
where now U{t) is the solution of the differential equation: 

jlJ{t) ^ -k{t)U{t) (12) 

with initial condition U{Q) — 1 and 

K(t) = - J dxJ{x,t)p{x,t) (13) 
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Following the same steps which led to Eq. (^) it is easy to show that the correlation functions p(xi,ti; ...Xn,tn) 
generated by S[J] 



1 ,5"S[J] 

p(Xi,ti] ...Xn,tn) = ^ 



(14) 

j=o 



E[J] 6J{xi,ti)...SJ{Xn,tn) 

reduce to those defined in (||). Analogously, it is useful to introduce the connected correlation functions defined by: 

(5"lnS[J] 



5J{xi,ti)...5J{Xn,tn) 



(15) 

J=0 



In the following, we will take the perturbative expansion (|^) as a definition of the fully interacting partition function 
where the reference n-particle correlation functions are written as functional derivatives via Eqs. ( p^ , p^ . Now the 
quantum to classical correspondence is complete: Having defined the (quantum) reference system by a hamiltonian 
Hr^ we compute the corresponding functional Sfl[J] via Eq. (pH). If we interpret this quantity as the grand partition 
function of a fictitious classical reference system in c? + 1 dimensions in an external field J(a;, t) coupled to a scalar 
order parameter ip(x,t) and at an effective temperature Td = 1, the exact partition function of the quantum system 
formally coincides with that of such a classical system with the additional non local coupling in the field ip{x,t) 
at two different space-time points. 

For future reference, let us briefly discuss the general relationship between the previously introduced imaginary 
time two point function ([l^) and the physically relevant dynamic correlations defined by: 

D{xi-X2,h-t2) = {e-^^e*^*iA/5(a;i)e^*^*ie'^*^Ap(x2)e-'^*^} (16) 

where Ap(x) = p{^)~ < pi^) Taking the Fourier transform in the space and time variables and inserting a 
complete set of basis states we get the usual Lehmann representation for the spectral function: 

lmD{k,uj) = — 51 e"''^'- I < s\pk\s' > \^S{uj-Es' + Es) (17) 

s,s' 

which obeys the detailed balance relation lmD{k,uj) = Iml?(— /c, — w) exp(/3w). The analogous expression for the 
imaginary time correlation function reads: 

F{k, c.„) = ^ ^ — \<s\pk\s'>\^ (18) 

s.s' 

where we took advantage of the periodicity of the two particle correlation function (^ in the imaginary time variable 
t = ti — t2 with period /3. As a consequence, the frequency variable w„ is now an integer multiple of 2tt/ fi. By 
comparing Eqs. ( p^ and ( |l^ ) we get the final relation, valid for systems with inversion symmetry: 



— lmD{k,uj) .'^ „ forcc;„^0 (19) 

-OO 1" ^ + 

This equation will be useful in order to relate the correlation function F{k,uj) to physical observables. 

In order to derive a set of renormalization gcoup equations for the model hamiltonian (|^), we follow the same 
strategy successfully adopted in classical systemaiJ: First we obtain a formal expansion of the Helmholtz free energy 
of the quantum model in powers of the two body interaction w{x). Then we define a sequence of intermediate systems 
( Q-systems) in which fluctuations characterized by wavevectors k < Q are inhibited. Finally we write a differential 
equation for the "evolution" of the free energy when the infrared cut-off Q is changed by an infinitesimal amount. 
Analogously to the classical case, the resulting equation is exact but cannot be written in closed form because it 
involves the knowledge of the two particle correlation function of the Q-systems. In fact, this is only the first equation 
of an infinite hierarchy of differential equations for the free energy and n-particle correlation function of the model, 
which form the Quantum Hierarchical Reference Theory of Fluids (QHRT). P-. 

Let us briefly sketch the basic steps necessary to derive the first QHRT equation: Details can be found in Ref.Ej 
via the quantum to classical correspondence previously discussed. The perturbative expansion of the logarithm of 
the grand partition function is directly obtained from Eq. (H) in terms of the connected correlation functions (|l5| ) of 
the reference system, denoted by F^. Next we perform a Legendre transformation which defines the Helmholtz free 
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energy A[p] as a functional of the order parameter profile p{x) = (3 ^ ((51n^/5J(x)) corresponding to a static external 
field J{x): 

A[p\ = -/3-MnS[J] + j dxJ{x)p{x) (20) 

The perturbation series of the free energy in powers of the interaction is then given by a diagrammatic expansion 
involving the pair potential (0) and the reference correlation functions (^5|). This series can be conveniently ordered 
according to the number of loops: the terms up to the one loop level reproduce the well known random phase 
approximation (RPA) to the free energy of the interacting system. In the homogeneous limit p{x) = p, the expansion 
reads: 

A = An + \vp^w{k = 0) - \vpw{r ^0) + V^Y.J ^ In [1 + Ffl(fc, LOn)w{k)] + £2 (21) 

Here V is the volume of the system, £2 represents the sum of all the diagrams with at least two loops, whose structure 
we do not need to detail, k is the wavevector and w„ — 27rn//3 are the usual Matsubara frequencies coming from the 
Fourier transform in the imaginary time direction, while w{k) is the Fourier transform of the two body interaction 
w{x). 

Now we define a sequence of Q-systems interpolating between the reference system and the fully interacting one 
by introducing a cut-off Q in the wavevector integration of every loop of the expansion ( ^l| ) . Clearly, for Q — > we 
recover the full expansion (^ij) to all orders while for Q ~* 00 all contributions to the free energy with at least one 
loop are suppressed and we get the mean field expression for the Helmholtz free energy: 

Aq=oo = Ah + ^Vp^w{k = 0) - ^Vpwir = 0) (22) 

Analogously, the two point function F{k,uJn) admits a perturbative expansion in powers of w(r) which can be again 
ordered according to the number of loops. Only the zero loop (chain) terms survive in the Q — > 00 limit reproducing 
the well known RPA result: 

1 + FR(fc, w„)w(fc) 

These results show that, through the sequence of Q-systems, thermal and quantum fiuctuations are gradually turned 
on and consequently, the free energy of the Q-system continuously connects the mean field approximation ( [2^j2^ ) to 
the exact expression corresponding to the fully interacting model. 

As already noticed for the classical casellJ, it is possible to write an exact differential equation describing the Q- 
evolution of the free energy Aq. The key observation is that in all the diagrams of the expansion (|2l]), including the 
terms collectively named £2, every loop contains at least an interaction bond. Therefore, limiting each momentum 
integration to the domain k > Q can be equally achieved by cutting off the Fourier components of the interaction w{k) 
at small momenta. Then, the Q-system can be equivalently defined as the system interacting via the pair potential: 

^^('^^^lo iffc<Q (24) 

The only caveat concerns the zero loop diagrams in the diagrammatic expansion of the free energy, which only depend 
on the zero momentum component of the perturbation potential WQ{k). According to the definition (p^, WQ{k = 0) 
identically vanishes for every Q 0: The zero loop contributions are then absent from the free energy corresponding 
to the potential (^), while were included in the previous definition of Q-system. Therefore, in order to get the free 
energy of the Q-system from the free energy of the hypothetical model interacting through the potential wq , wc must 
explicitly add the zero loop contributions, which have been already evaluated in Eq. (p^). Similar considerations hold 
for the two point function of the Q-system Fqlk) which differs from that corresponding to the interaction wq because 
of the simply connected (chain) diagrams. However, these terms are proportional to the A;*'' Fourier component of 
the interaction 'WQ{k) which, according to the definition (p4[), coincides with w{k) for every k > Q while vanishes if 
k < Q. Therefore, the chain contribution needs to be added to the two particle correlation function of the system 
with interaction wq only for k < Q. As a consequence, the two point function of the model with interaction wq, 
which we denote by FQ{k,uJn) is discontinuous in k space on the k = Q surface and can be expressed in terms of the 
two point function of the Q-system FQ{k) as 
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FQ{k,ujn) if fc > Q 

fQ(fc,^„)^<! fQ(fc,^„) .^^^^ (25) 

1 - FQ{k,ujn)w(k) 

This one to one correspondence between Q-systems and models interacting via the potential ( p^ is extremely useful 
in order to obtain the QHRT evolution equations. According to the previous discussion, the change in the free 
energy of the Q-system corresponding to a change in the cut-off Q ^ Q — 6Q can be evaluated by means of the 
perturbative expansion ( pl| ) in which the "reference" free energy is now identified as Aq and the perturbation potential 
as [wQ-sQ{r) — 'WQ{r)]. Such a potential has Fourier components only in the infinitesimal momentum shell defined 
hy Q — SQ < k < Q. In the limit SQ 0, by keeping only terms 0{6Q) we just select the one loop diagrams which 
have been explicitly summed up in Eq. (pi|), while all the many loop contributions £2 are at least of order (SQ)"^, 
i.e. rigorously negligible. The zero loop terms are trivial and have been automatically included form the beginning 
through the initial condition (p2|). The resulting differential equation for the evolution of the free energy density is 
then: 



dQ\V J 2/3 



E/ -^,H^-FQ{k,u.^)w{k)] (26) 



where the momentum integral is limited on the shell k — Q and then reduces to the integration over the solid angle 
rik- 



This evolution equation is exact and must be supplemented by the initial condition (22) at Q = 00. However, the 
presence of the two point function FQ{k, uin) at right hand side prevents the possibility to obtain the thermodynamics 
of the model without introducing some sort of approximation. In the following we will discuss a closure to Eq. ( [2^ ) 
based on a parametrized form of the two point function which is related by an exact sum rule to the Helmholtz free 
energy Aq of the Q-system. The justification of this closure relation requires, however, an analysis of the physical 
properties of the model we want to investigate and therefore we first have to specialize this rather general presentation 
of the QHRT equations to the case of the antiferromagnetic Heisenberg model which will be thoroughly investigated 
in the following Sections. 



III. QHRT EQUATION FOR THE ISOTROPIC HEISENBERG MODEL 



In this Section, the QHRT formalism just developed will be applied to the study of the isotropic antiferromagnetic 
Heisenberg model on a iV-site cubic lattice in arbitrary dimension d and spin S. This is an extremely interesting 
many body model which allows to study the interplay of thermal and quantum fluctuations in determining the 
phase diagram of the system and the effects of spatial dimensionality on the occurrence of long range orderllll. In 
the 5' — *■ 00 limit, quantum fluctuatioiis are suppressed and tiie system behaves classically. Thija regime has been 
extensively studied by aeries expansionsEI, numerical simulationsB, renormalization group techniquesB and approximate 
theoretical approaches^. Much less is known at finite values of the local spin S: accurate simulations have become 
available only recently while series expansions are known to lose accuracy both at low temperature and in the critical 
region. 

Since we want to investigate the paramagnetic-antiferromagnetic transition, it is convenient to consider the presence 
of an external magnetic field h, parallel to the z axis, which couples to the order parameter, i.e. to the staggered 
magnetization: 

H = J Sr-Sr, (27) 

<R,R'> R 



Here J > 0, g is the antiferromagnetic wavevector of components gi — tt, while the spin operators S satisfy the usual 
commutation relations 



(28) 



The lower and upper indices label respectively the site on the lattice and the component of the spin, while eijk is 
the totally antisymmetric tensor. With respect to the class of hamiltonians examined in the previous Section, the 
Heisenberg model requires some straightforward extension because i) the system is defined on a lattice, which is not 
isotropic; ii) the order parameter is now the staggered magnetization: a three component vector. As a consequence. 
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the correlation functions acquire spin component indices and Hi) the presence of a staggered magnetic field breaks 
translation invariance by one lattice site. 

The reference system is chosen to be defined by the non interacting hamiltonian 

HR = -hJ2e'^'^SI^ (29) 

R 

whose properties can be easily obtained. For instance, the grand potential as a function of the number of sites N, 
inverse temperature /3 and staggered field h reads: 

log Er = Nln [smh{f3hS) coth(/3/i/2) + cosh{(3hS)] (30) 

The staggered magnetization m can be found by differentiating ( ^ ) with respect to l3h and the Helmholtz free energy 
then follows by Legendre transform (pO[). However, the explicit inversion formulae can be found analytically only in 
the 5 = 1/2 case where 

to) 

(31) 

In order to derive the renormalization group equation appropriate for the magnetic transition in this model, we now 
have to suitably define the sequence of Q-systems whose role is to gradually introduce fluctuations into the system. 
The critical fluctuations occur at a wavevector k = g and then it is natural to define a sequence of systems in which 
the Fourier components of the potential ^(k) 

d 

w{k) = 2 J7(k) 7(k) cos ki (32) 

i=l 

are suppressed in a neighborhood of k = g. A convenient choice is to limit the domain of integration in k space by 

|7(k)| < Vd^~Q^ (33) 

with Q G (0, d). As Q — > all Fourier components are included and we recover the fully interacting system. At Q = d 
the free energy density a ^ A/N includes only the reference term and the zero loop contribution: 

aQ=d = CLR- dJrn^ (34) 

The two particle correlation function is represented by a a 3 x 3 matrix which, due to the spin symmetries in the 
Hamiltonian under rotations around the z axis, has the following general structure in Fourier space: 



N 



{\ + m) ln(i + to) + (i - to) ln(i 



(3h = In 



+ TO 



F 



pxx pxy Q 

pyx pvv 
F^^ 



(35) 



The three independent entries are given by: 

i^""(ki,c^i;k2,c^2) =-F^^(ki,a;i ;k2, 6^-2 ) =F""(ki,wi)/3iV(5(wi +6^-2) 5(ki +k2) 

F-y{\<ii,LUiM.^2) = -Fy%lii,Lo^-M.uJ2) - F^nki,a;i)/3^'5(a;i +cc>2)^(ki +k2 -g) 

F^' (ki , wi ; k2 , 6^2 ) = F^' (ki , wi ) /3 N5{ux + wa) '5(ki + kj ) (36) 

where k belongs to the first Brillouin zone and uo is an integer multiple of 27r//3. Because of the presence of the external 
field, the translational invariance is broken and the two point function is not proportional to (5(ki +k2). Nevertheless, 
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due to the symmetry of the hamiltonian under a translation of one lattice site and a simultaneous rotation of tt 
around x axis, the diagonal elements and the off-diagonal elements are respectively proportional to (5(ki + and to 
5(ki + k2 — g) as stated in Eq. (|36|). At Q — d, the summation of the zero loop diagrams explicitly gives: 

<3='^^ ' ' m-^Lj^ + Mi - ^(k)2 

F^L,(k,c.) = ^^ (37) 
^ ^11 + iy(k) 

where fi± and are functions of temperature and staggered magnetization m given, for S = 1/2, by: 

Pfi±_ = 2m-^ tanh"\2m) /3^|| = 4(1 - W)"^ (38) 

Having determined the mean field form (i.e the zero loop contribution) of the relevant quantities, we now turn to 
the inclusion of fluctuations via the definition of the QHRT renormalization group equation. Following the procedure 
outlined in the previous Section and taking into account the definition ( |33| ) of Q-systems, we find the exact evolution 
equation for the free energy density aq = Aq /N of the antifferromagnetic Heisenberg model which should be solved 
with the appropriate initial condition (|3^) at Q = d: 



da, 



Q 



dQ hp^j ^Wd^ - 7(k)2 - Q) { ln(l - ^^^^(k, c.„)«)(k)) + (39) 



In 



(1 - ^^^-(k, c.„)iZ,(k))(l + F^-(k + g, c.„)il;(k)) - F^^(k, Un)F^J{\^ + g, u^)w{\^f] } 



The summation of the one loop (ring) diagrams, necessary to obtain Eq (p9|), must be carried out by carefully exploiting 
the frequency and momentum conservation rules (|3q). The delta function present in the d-dimensional integral (pS) 



limits the integration domain to the surface Sg defined by 7(k) = cP — Q^. This equation, although formally 
exact, is not closed because the evolution of the free energy depends on the unknown magnetic structure factors of 
the model Fq (k, oj) at a generic value of the cut-off Q while these quantities are known only at the beginning of the 

integration, i.e. at Q = d where they have the form shown in Eq. (|3^). If we simply assume that (k, w) is not 
affected by fluctuations, the QHRT equation reproduces the the standard RPA approximation. A better treatment of 
fluctuations needs a more flexible parameterization of the correlation functions. This delicate point requires a detailed 
analysis of the properties of the (imaginary time) dynamic structure factors of the Heisenberg model, which will be 
carried out in the next Section. 



IV. TWO BODY CORRELATIONS 



In this Section we exploit some known feature of quantum antiferromagnets in order to derive constraints on the 



formal structure of the two point functions i^g^ (k,cj„) required to close the QHRT evolution equation (|39|). 

As a first step, let us recall the relationship between the dynamic correlations of the Heisenberg model, defined in 

physical frequency, and the imaginary time functions which enter the QHRT formalism. The general correspondence 

has been already obtained in Eq. (|l9|): Here we point out that the analytical form of the mean field correlation 

functions (p7|), which serves as initial condition, reproduces the single mode approximation in the dynamic structure 

factor: 

ImS'^^(k,t^) 



1 — exp(— /i_L + w(k) 



Im^-(k,.) = ^^(^^ (40) 
p ^11 + w(k) 



where 



tk = m. hil - w(k)2 (41) 
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can be recognized as the dispersion relation for spin waves in this approximation. The elastic sum rule immediately 
gives the corresponding form of the equal time structure factor via integration over the frequency axis: 



Efc 1 



' ' ' 2tanh(/3efc/2) + w(k) 

^ii(i^) = ^E^^^(k--0 4^^r;:^ (42) 

These expressions show that both the dynamic and static structure factors ( ^ , ^2| ) coincide witk-Lhose usually adopted 
in the interpretation of neutron scattering data, in the limit of negligible spin wave dampingl3. This suggests that 
the simple analytical structure of Eqs. l pT\ ) qualifies as a good starting point for obtaining a closure relation to the 
QHRT equation. However, the parameters fi± and appearing in the correlation functions (|37| ) will be modified by 
fluctuations and cannot be identified with their mean field estimate ( |38| ) down to Q = 0. These quantities can be 
directly related to important physical properties of the model: The definition ( |l5|) of n-point function shows that the 
zero frequency value of i^"(k, a;„) coincides with the uniform (staggered) susceptibility x° ix") when k = (k = g) 
and therefore 

xl = F-^0,0)^if,^ + 2Jd)~' 

X;i=i^^nO,0) = (Ai||+2Jrf)-i (43) 



^S = (£)" = (0)" = ^"fe''')-(''ii-"'"- <«) 

Here, we have conveniently expressed the staggered susceptibilities as thermodynamic derivatives of the external field 
h and of the free energy density a with respect to the staggered magnetization to, using the isotropy of the model in 
the spin variables. At the critical point and in the whole coexistence region, both the longitudinal and the transverse 
susceptibilities divergea, implying that = = 2Jd. This feature leads to a gapless spin wave spectrum, via Eq. 
(p]), with linear behavior near k = (and k = g): ~ c^fc. In the approximation (p7|), the associated spin velocity 
is proportional to the spontaneous staggered magnetization: Cs — Jm\/4d. At zero temperature, the static transverse 
structure factor (^) is linear in k around k = in the symmetry broken phase and the prefactor defines the spin 
stifltness pjE3: 

S^{k)^^k (45) 



Within a parameterization of the form (37), the stiffness can be expressed in terms of the spin velocity: ps = 
c'^/{4:Jd) — Jm?. Remarkably, the hydrodynamic relation x°l'^s = Ps, which is believed to hold at zero temperature, 
is satisfied even in this simple approximation. 

A known sum rule which is worth mentioning regards the on site value of the static structure factors, which is 
bound by the spin normalization condition: 



/3 



at any temperature. If we choose to determine the parameters /j. j and by use of the susceptibility sum rules (|44|), 
in close analogy with the procedure adopted for classical modelsEfl, the additional relation (^6|) cannot be generally 
satisfied by the simple form (|3^) . This deficiency is due to the presence of an incoherent part in the actual dynamic 
correlations which has been neglected in our single mode approximation. 

A useful check on the accuracy of the parameterization (p^) comes from the zero temperature limit, where quantita- 
tive approximations are available, like Spin Wave Theory (SWT). Although SWT is formally an expansion around the 
classical limit, even the simplest truncations in the parameter 1/ S are known to be quite accurate even for 5* = 1/20. 
To lowest order, the spontaneous staggered magnetization attains its maximum value m = S and the correlation 
functions in imaginary time can be easily calculated in vanishing external field h: Remarkably, the result has precisely 
the mean field structure (^) with fi± = = 2Jd, as expected in the symmetry broken state. However, it is known 
that quantum fluctuations affect the zero held magnetization, the spin wave spectrum e^-, the stiffness ps and the 
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uniform susceptibility x°\_ ■ These quantities are indeed modified to first order in swtEI, although the form of (Ul^ 
remains unaltered and only the prefactor is changed, leading to a simple renormalization of the spin wave velocity 





f d'^k 


[1 


- es«,(k)] 


'-hi 


f d'^k 


[1 


- eswOi)f 








''2SJ 


f d'^k 


1 




' {2Tr)d 







Cs ~ Zc JSv Ad ; Zc 
Ps = Zp JS^ ; Zp 

,47) 

where the momentum integrals are restricted to the lattice Brillouin zone and 



e.„(k) = ^l-d-2^2 (48) 

A better representation of the dynamic correlations which, keeping the single pole structure (|4^), is able to reproduce 
the first order SWT results @ can be obtained by including two more adjustable renormalization amplitudes ^ in 
the functions i^'-' (k, a;„): 



fi± - u;(k) 



m 



F-(k,c.) = C||^^ (49) 
Mil +w(k) 

At zero temperature, the full one loop SWT corrections ( ^7|) are reproduced by the form ( ^9| ) with the choice (± = Z^ 
showing that such an ansatz for the dynamic correlation functions may improve the simple mean field description 
of Heisenberg antiferromagnets. However, in order to implement this generalization we need additional sum rules 
for determining self-consistently the parameters C± and besides the already mentioned Eqs. (^4|) for fj,± and fi\\. 



A natural choice would be the analogous Eqs. (43), which control the renormalization of the uniform susceptibility. 
However, we will not explicitly pursue this extension in the present work, but rather confine our analysis to the 
simplest choice Cj- = C|| li i-c keeping the same structure of the dynamic correlations at the mean field level 
(p7|), just allowing for "mass renormalization", according to Eqs. (^). The single mode structure of the dynamic 
correlations, besides reproducing the lowest order spin wave result at X = 0, has also been proven quite accurate in 
numerical studies on the Heisenberg model both via series expansionsEJ and Monte Carlo simulationsu. 

As noticed by Hertz in his pioneering work on phase transitions in quantum systema23, a proper treatment of critical 
phenomena in quantum mechanics directly provides predictions on dynamic scaling in the system. The same feature is 
clearly present in our microscopic approach which requires the form of the dynamic structure factor. Unfortunately, in 
the approximation we are examining, the single pole form ( |4C| ) does not allow for anomalous scaling in the dispersion 
relation (w ^ k^) at the critical point. The functional dependence on k and lo in Eq. ( p] ) correctly predicts a linear 
spin wave dispersion to ~ Csk in the symmetry broken state at all temperatures T < Tc- On approaching Tc, however, 
the spin wave spectrum vanishes as the order parameter ek ^ rn ^ where t — {Tc — T)/Tc and (3 is the appropriate 
critical exponent. Therefore, exactly at Tc, the spin wave spectrum fades away and above Tc only a delta contribution 
at w = survives in the paramagnetic dynamic correlations. Such a description of dynamic critical behavior provided 
by the approximate form ( ^ ) is clearly rather schematic and unsatisfactory. However, the way the spin velocity 
vanishes Cg ^ is consistent with the predictions of dynamic scaling Cg ~ if the dynamic critical exponent is 
given by z = d/2. In fact, by inserting the power law behavior of the correlation length ^ ~ t^'^ near the critical 
point, we obtain Cs ^ iv{d-2)/2^ which gives Cs ~ via standard hyperscaling relations. The latter equality holds 
because in our approach the correlation exponent = Q due to the analytic k dependence of the correlation functions. 
Remarkably, the result z = d/2 agrees with the general expectations of dynamic scaling for antiferromagnets: The 
appropriate equation of motion corresponds to "model G" in which the order parametar does not commute with the 
hamiltonian but a conserved quantity (i.e. the uniform magnetization) is also presentEj. 

The discussion presented in this Section supports the choice of a simple single mode form as a reasonable approxi- 
mation of the actual dynamic structure factor in Heisenberg antiferromagnets and leads us to propose Eqs. ( p7| , p4[ ) as 
a closure to |tlie exact renormalization group equation (^, closely following the strategy successfully adopted in clas- 
sical modelstil. In fact, the thermodynamic relationship (E3) between correlations and derivatives of the free energy 
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holds independently of the form of the interactfon ty(k) and then it is valid for every Q-system, which corresponds to a 
potential WQ(k) characterized by a cut-off Q on its Fourier components. Therefore, the dynamic correlation functions 
of the Q-systems Fq can be parametrized by the same form ( |37| ) used at Q — d with Q-dependent parameters /i_LQ 
and /i||Q defined by the analog of Eqs. (44): 



m 



dm I 



dar 



dm? 



= i^Q^(g,0) = (/X||Q-2Jd)-i 



(50) 



When the approximate form (^7|, pO| ) of the dynamic correlations is inserted in the exact evolution equation (|3^) we 
get a closed equation for the evolution of the free energy of the Heisenberg model. The frequency sums in Eq. ( |39| ) 
can be analytically evaluated by Euler summation formula while the momentum integration can be trivially carried 
out with the result: 



daq _ Dd{Q) 



4 In 



sinh (i/Jm^j^g) 



2/3 1 " [sinh (^^Pm^nlQ - AJ^{d^ - 
where the "density of states" Dd{Q) is defined as 



In 



^IIQ 



Aif|Q-4j2(d2_Q2) 



(51) 



(52) 



This evolution equation will be analyzed in the following Sections. 

Finally, it is interesting to consider the limiting form of our equation ( pT| ) when S* — > cxd, i.e. when the quantum 
Heisenberg model reduces to a classical 0(3) vector model. Clearly, we have to rescale all the dimensional quantities 
by suitable powers of S in order to get a finite result. In particular, the magnetization m — s- Sm while energy and 
temperature must be scaled as a ^ S^a and (3 S^'^p. As a consequence, equation (51) simplifies and we get 



daq 
~dQ 



Dd{Q) 
2/3 



2 In 



MiQ-4j2(d2_Q2 



+ ln 



/x2Q-4J2(d2_Q2) 



(53) 



which in fact coincides with the equation obtained in the framework of HRT for a classical vector mode 
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V. LONG WAVELENGTH ANALYSIS 

Let us study the behavior of the QHRT evolution equation (|5^) in the critical region. This investigation is par- 
ticularly instructive because it shows how universality and renormalization group (RG) come about in the QHRT 
framework. We will consider the critical region of the paramagnetic-antiferromagnetic transition, where the physical 
staggered susceptibilities (xj_, xf|) diverge. First we remark that our ansatz for the two point function (^) is analytic 
in k. As a consequence, the correlation exponent predicted by our approach is 77 = in any dimension. It is known 
that such a result is correct only above the upper critical dimension d = 4. However, due to the smallness of 77 in 
d = 3 (ry ~ 0.04u), we consider such an approximation acceptable also in the physically relevant applications. 

In order to study the critical behavior predicted by our approach, we now focus our attention on the last stages of 
the evolution {Q ~ 0), i.e. on the long wavelength form of equation (^l|). In fact, we know that the order parameter 
fluctuations of wavevector k close to g are inhibited by the presence of a cut-off Q and the free energy of a Q-system 
cannot develop singularities at finite Q, but only in the Q limit. Eqs. ( pO| ) show that, in this regime, both /ij^Q 
and /i||Q tend to 2dJ . Moreover, being the critical point at m = 0, our analysis will be limited to the to ^ region. 
The asymptotic form of the QHRT equation (^ij) is then: 



daq 
'dQ 



2P{2TTVdy 



Q 



d-l 



2 In 



dr 



In 



J 



(54) 



where fi^ is the surface of the unit c?-dimensional hypersphere. Remarkably, this limiting form of the evolution 
equation coincides with the one appropriate for a classical 0(3) vector model in d dimension, i.e. with the long 
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wavelength limit of Eq. (|5^). This, confirms that Quantum Mechanics becomes irrelevant in the neighborhood of a 
(finite temperature) critical pointer. 

In order to exploit the renormalization group structure hidden in Eq. (p4h we suitably rescale the variables: 



H 



2(3J{2TTVdY_ 



'1/2 



Q 



-{d-2)/2. 



2/3(27r^/rf)'i 



Q [aqim) - aqim = 0)] 



In this way, the explicit dependence of the equation on Q is eliminated and the evolution equation becomes: 



H ^ dH + z'^—J^H' 



In 



H" 



1 



H"{Q) + 1 



2 In 



H'/z + 1 



iJ"(0) + 1 



(55) 



(56) 



where dots represent derivatives with respect to I = \n{d/Q) and primes derivatives in z. This-equation is formally 
identical to that obtained by NicoU et al. in an approximate renormalization group formulationES and, as such, it has 
already been studied by means of standard RG techniques. .-This analysis showed that the critical point is present 
only for d > 2, in agreement with Mermin- Wagner theoremllZl. Note that precisely at T = 0, the QHRT equation 
( pl| ) does not reduce to the asymptotic form ( p4|) and therefore quantum critical points would require a diffeijent 
long wavelength analysis. Eq. (|5^) properly embodies critical scaling and hyperscaling and, as shown in Ref.E3, it 
correctly gives non classical critical exponents to first order in the e = A ~ d expansion. For instance, the staggered 
susceptibility exponent at m = (x|*| ~ t~^) is; 
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1 



22 



0(e^ 



(57) 



The value of 7 predicted by our approach in d = 3 can be found by integrating the fixed point equation associated 
to ( |56| ) and solving the corresponding eigenvalue problem. The nuiaedcal result gives 7 = 1.65, to be compared with 
the accepted value for a three component order parameter 7 = 1.3900. The other critical exponents can be obtained 
from 7 by scaling, recalling that = in our approximation. 

Interestingly, a similar RG analysis can be carried out on the same asymptotic QHRT equation (^4|) in order to 
study the divergence of the susceptibility on approaching the coexistence curve at T < Tc- By rescaling m in the 
neighborhood of the spontaneous magnetization TOx as 



d flrl 



H 



-i -1 



_2/3{2TTVdY 

the evolution equations at long wavelengths reads: 

H = dH + z{2- d)H' 



Q [aQim) ~ aQ{my_)] 



In 



H" 



H"{0) 



2 In. 



H' + 1 
H'{Q) + 1 



(58) 



(59) 



At coexistence, the free energy approaches a fixed point which corresponds to a stationary solution of of Eq. (|5S| ), 
regular in the whole interval z S (0, cxd) and satisfying the boundary conditions H{0) = and H{z) z'^/(''~^'as 
z — > +00. The latter requirement is necessary in order to match with the solution outside the coexistence region (i.e. 
at m > mx). As a consequence, the staggered susceptibilities display a non-analytic behavior when d < 4: 



X± 



(m-TOx)"'/*'"'^ 



(m-mx)-^'-')/^'-') 



(60) 



in agreement with the predictions of a spin wave analysis in the classical 0(3) modeld. 

Let us conclude this section showing how it is possible to recover, within the QHRT formalism, also the RG equations 
for the coupling constants of the quantum non linear sigma model [QN LaM). The long- wavelength behavior of the 
d dimensional Heisenberg antiferromagnet in zero external field-^ commonly believed to be described by a d -I- 1 
dimensional QNLcrM, at least in the low temperature regimecZl. This is an effective field theory in which the 
additional temporal variable, defined only in the interval (0,/3), reflects the presence of quantum fluctuations. The 
QNLaM action is conventionally written in terms of a field S},{r,t) as: 
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(61) 

where the local constraint |n| = 1 is understood and Ps and Cg are the bare (i.e. mean field) spin stiffness and velocity. 
The phase diagram of the model depends on temperature and on the dimcnsionless coupling constant: 

9 - (62) 

where A is the ultraviolet momentum cut-ofF, which is of the order of the reciprocal of the lattice spacing. The 
QNLaM , as a low energy iheory for Heisenberg antiferromagnets, has been extensively studied by Chakravarty 
Halperin and Nelson (CHN)El by weak coupling RG analysis. The one loop equations describing the evolution of the 
couplings constants are: 



^ = {l-d)g+^g'coth{g/2t) 



dl\t) 



-f (63) 



Here, is the length rescaling factor and Kd is a geometrical coefficient (i^J^^ = 2'^~^7r'^/^r(c?/2)). The equations are 
written for the running coupling constant g{l) and for the dimcnsionless temperature scale t{l) whose initial condition 
\s t ~ A'^^^/(/3ps)- For d < 2 there are no finite temperature fixed points, while at T = and for c? > 1 there is a 
non trivial fixed point given by gc — 2/ Kd{d — 1) describing a quantum transition. Depending on the initial value of 
the coupling constant the RG flow drives g either to zero or to infinity. In the former case the physics is governed 
by the low temperature fixed point and then the system is ordered. Assuming that the two dimensional Heisenberg 
model has a symmetry breaking ground state, from the solution of the RG equations (^3|), the following behavior of 
the correlation length is obtaineda: 

^{T) ^ eoe^ (64) 

where ps = Psidc — d)/ 9c is the renormalized spin-stiffness at T = 0. 

The weak coupling RG analysis predicts several regimes where quantum and thermal fluctuations play different 
roles. However, quantitative estimates based on low energy effective theories are usually rather inaccurate because 
short wavelength fluctuations are neglected in these approaches. Here we want to show that the QHRT equation, 
reduces, at low energy, precisely to the one loop RG equations studied by CHN (^). However, QHRT provides a 
description of spin fluctuations over all lengthscales maintaining the advantages of a fully microscopic theory. 

In the low temperature region, longitudinal fluctuations can be neglected because they are characterized by a less 
s ing ular behavior and therefore, in this regime, the last term in Eq. (pTl) may be dropped. If we now differentiate Eq. 
( pl[ ) with respect to the magnetization, we get an equation for the evolution of the staggered fleld hQ{m) at fixed m: 



dh^ 



Q 



dQ 



coth (^m/ixQ 1 



Dd{Q) 

cothf^^MiQ-4P(d2-Q2) 



2 V ^ } V^Q-4j2(d2_Q2) 



(65) 



This equation implicitly describes how the spontaneous magnetization is reduced from its mean field estimate when 
fluctuation are included. In order to exploit this information, it is useful to define a (^'dependent spontaneous 
magnetization mg by requiring that /iq(toq) = 0. Such a quantity satisfies an evolution equation which can be easily 
deduced from Eq. (|65|). At long wavelengths (i.e. keeping the leading terms as Q ^ 0) the Q dependent staggered 
magnetization in zero field obeys a closed equation: 



-Kd\^^J coth(/3QmQ) (66) 

In order to make contact with the RG picture based on the QNLaM effective action, we now have to relate the 
QHRT variable mg with the running coupling constant g{l) appearing in Eqs. (^^. The correspondence between 
the length rescaling factors is clearly I — \n{d/Q), while the definition of the couphng constant g in terms of physical 
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quantities, like stiffness and spin velocity, given in Eq. (|62[) suggests the proportionality between and the 
spontaneous magnetization. In fact, as discussed in Section (|V|), QHRT provides explicit relations between ps, Cg 
and the magnetization at coexistence. It is therefore rather natural to introduce an "effective coupling constant" 
gg = V4d{Q/\/d)'^~^mQ^ which, when substituted in Eq. (|66|), satisfies an equation formally identical (to order g^) 
with the standard one loop RG form (^3|). As an example, we plot in Figure the RG flux of gq as obtained by 
the integration of the full QHRT equation in d = 2 ( pT| ) via the previous definition of gq . We clearly see the effect of 
the unstable zero temperature weak coupling fixed point while, for the nearest neighbor Heisenberg model, the other 
fixed point (gc — 47r), governing the Quantum Critical regime, has no effect on the RG trajectories. This suggests 
that the Heisenberg model always remains in the Renormalized Classical regime in agreement with an independent 
analysis based on the pure quantum self consistent harmonic approximation (PQSCHA)E3. 



VI. NUMERICAL RESULTS 



Equation (|5l|), together with the thermodynamically consistent closure relations (Q) and the appropriate initial 
condition at Q = d (^) has been solved numerically in two and_three dimensions and several values of the spin S. 
Following the analogous treatment of HRT in classical systemalj we fcst wrote the partial differential equation in 
quasilinear form and then used a fully implicit finite difference schemata. A careful implementation of the numerical 
algorithm is necessary, because the solution to the QHRT equation develops singularities as Q ^ 0, when a phase 
boundary is crossed. This is a consequence of the special treatment of long wavelength fluctuations provided by 
QHRT which is able to reproduce rigorously flat isotherms inside the coexistence region, as already shown in classical 
systemsca. Our computations were carried out with a mesh of few hundred points along the magnetization axis, 
which already allows for a good accuracy. The high temperature asymptotic behavior of thermodynamic quantities is 
correctly reproduced by QHRT in arbitrary dimension d and for every spin S, as can be easily shown by analysis of 
the evolution equation (pl|). Before discussing the numerical results in three dimensions we briefly comment on the 
calculations performed on the square lattice. In this Section we set J = 1. 



A. The two dimensional case 



As already noticed, our equation correctly predicts the absence of spontaneous symmetry breaking at finite tem- 
perature for every choice of the spin S. However, the effective coupling constant fiow shown in Fig. illustrates 
the strong attractive nature of the unstable zero temperature fixed point and allows to give a precise estimate of the 
ground state spontaneous magnetization. For S ~ 1/2 we obtain mx ^ 0.35, rather close to the value m ~ 0.307 
obtained by numerical simulationaD. Weak coupling RG analysis predicts that the correlation length at m = diverges 
exponentially on approaching zero temperature with a coefficient related to the spin stiffness (|^) . The same behavior 
is found by our numerical results as shown in Fig. j2j. This result is not unexpected because we have already shown 
that QHRT contains all the information provided by the one loop RG approach. The estimate of the spin stiffness 
coming fram a linear fit of the data in Fig. j2| gives ps ~ 0.12 for S = 1/2, to be compared with the accepted value 
Ps ^ 0.180. Our result is in fact consistent with the relationship ps = discussed in Section (IV) which follows 
from the adopted closure. Therefore we are led to conclude that the discrepancy between the QHRT and Monte Carlo 
(MC) in the estimated values of ps is probably due to the crude approximation for the two point functions we have 
adopted in this calculation, as anticipated in Section (^V|), see also Eq. (p9|). 

The zero field specific heat Ch as a function of temperature is shown in Fig. |^ together with a comparison with 
MC data for S — 1/2lJ. The non monotonic behavior, even in the absence of a phase transition, is just a consequence 
of the asymptotic decay of Ch at high temperature {Ch oc T~^) and to the vanishing of Ch at T = 0, as required by 
thermodynamics. 



B. The Heisenberg model on a cubic lattice 



In d = 3, spontaneous symmetry breaking at finite temperature does occur in the Heisenberg model. As a 
consequence, thermodynamics forces the Helmholtz free energy to be rigorously flat inside the coexistence region. 
This feature is usually lost in the approximate treatments of the model and may be recovered a posteriori by Maxwell 
double tangent construction. As already noticed in the application of HRT to classical statistical models, our approach 
instead is able to predict rigorously flatisotherms in a region of the phase diagram, thereby allowing for a unambiguous 
determination of the coexistence curvecj. As an example, in Fig. B we show the behavior of the equation of state h{m) 
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at two temperatures: above and below for several values of the evolution parameter Q. We see that, although the 
initial condition at Q = 3 is in both cases below the mean field critical temperature (i.e. corresponding to a negative 
mean field transverse susceptibility), Xa} — h/m is strongly renormalized by fluctuations and for T > Tc it becomes 
positive and saturates as Q — *■ 0, while for T < eventually sticks to zero. Therefore, during the evolution, the 
region characterized by negative susceptibility is gradually removed (for T > Tc) or replaced by a flat isotherm (for 
T < Tc). ThiaJpehavior can be analytically understood from a study of the asymptotic equation, in close analogy to 
the Ising casec3. 

The full|-coexistence curve for S —rl/2 and S* = cx) is reported in Fig. Q together with the results of series 
expansionsa and numerical simulationsa. The zero temperature limit of the spontaneous magnetizalion gives my, = 
0.43 in our approximation, to be compared with the value 0.42 from second order spin-waves theorycl. The predicted 
critical temperature is shown in Table I for several values of the spin S. When available, also estimates from Monte 
Carlo, series expansions or cumulant expansions are given. In Fig. |^ the speciflc heat as a function of the temperature 
is plotted for S = 1/2 and S = oo together with classical MC simulations and recent QMC datacl. Note that QHRT 
correctly reproduces the different zero temperature limit of the quantum [Cr 0) and classical [Ch NJ) systems. 
Again, the agreement is quite satisfactory, being comparable with simulation uncertainties. We remark that QHRT 
correctly predicts a cusp- like behavior at Tc, although the corresponding critical exponent a < is underestimated 
by our approximation. These small differences, however, cannot be appreciated on the scale of the plot. Finally, Fig. 
W shows the change in the shape of the coexistence curve due to the effects of quantum fluctuations: going from 
the classical {S — oo) to the quantum S ~ 1/2) case, the magnetization curve rounds off and the linear behavior at 
low temperature disappears, in agreement with SWT which predicts a quadratic temperature dependence of mx as 
T ^Q. 

In summary, we believe that even within the crude simple pole approximation adopted in the present study, QHRT 
is able to get quantitative agreement with recent numerical data on the thermodynamics of quantum antiferromagnets 
in d = 3. The error can be estimated in few percents on the critical temperature, spontaneous magnetization and 
specific heat. Larger discrepancies can be seen in the shape of the coexistence curve: such a disagreement is rather 
unexpected and more accurate simulation data are needed before reaching a definite conclusion on this delicate issue. 



VII. CONCLUSIONS 



In this paper we have introduced a new theoretical tool for the study of phase-transitions in quantum systems. Such 
an approach is based on the use of the Hierarchical Reference Theory (HRT)Ej, a theory developed in the context 
of classical fluids able to correctly describe the system in a wide region of the phase diagram, both near and far 
from the critical point. Fluctuations over different lengthscales are gradually included thereby reproducing the full 
structure of Wilson renormalization group theory in the critical region. The possibility of applying HRT also in the 
context of quantum systems (QHRT) follows from the observation that the grand partition function of a wide class of 
quantum models can be formally mapped onto a classical partition function (in d-\- 1 dimensions) to which HRT can be 
directly applied. The advantage of our formulation is that since this mapping is exact, the microscopic informations 
are correctly included into this formalism and only in the vicinity of the critical point and at long wavelengths the 
peculiar features of the model disappear and the evolution equations acquire a universal form. QHRT is a general 
approach which can in principle be applied to a large class of systems including fermions, bosons and spins. 

The validity of this theory has been studied in the specific case of the isotropic Heisenberg antiferromagnet. This 
choice is dictated lw,the large amount of available results from other theories and numerical methods, which include 
mean field theorissp3, renormalization groupO and path integral approachcsli3, series expansionaJ or Quantum Monte 
Carlo simulationsa. In order to obtain a closed and tractable equation out of the QHRT formalism, we analyzed a 
simple approximation which amounts to parameterize the dynamic correlations in a single pole form. Even within 
this approximation scaling and hyperscaling are satisfied close to the critical point, the critical exponents are correct 
near four dimensions and, in d = 3, take non classical values slightly larger than the accepted ones. Remarkably, the 
weak coupling renormalization group equations for the Non Linear Sigma Model are recovered at low temperature 
even in d = 2, thereby reproducing the correct asymptotic scaling as T — > 0. 

The theory has been tested against numerical data on a simple cubic lattice with nearest neighbors interaction 
and for several values of the spin. The results are in close agreement with available simulations for several quantities 
like the location of the critical point, the spontaneous magnetization and the specific heat. In c? = 2 the comparison 
with simulations is less satisfactory even if QHRT is able to reproduce the qualitative behavior of the model. The 
discrepancies we found are related to the single mode approximation adopted in the closure of the equation and may 
be overcome by use of a more flexible form for the dynamic structure factors. 

In this paper our attention has been focused on the Heisenberg model since it allows for several interesting gen- 



15 



eralizations, like the introduction of anisotropies or competing interactions. In these cases, quantum and thermal 
fluctuations may play a more fundamental role in modifying the topology of the phase diagram. A detailed study of 
the Heisenberg model with easy plane anisotropy is under way: the results will shed light on the thermodynamics of 
the crossover between two and three component order parameter, with direct experimental implications. An extremely 
interesting application of the QHRT formalism to fermionic systems regards the Hubbard model, whose hamiltonian 
has in fact the general structure of Eq. (|l|) . The natural choice, as a reference system, is the free Fermi gas on the 
lattice in an external magnetic field, while the Hubbard interaction couples the on site electron densities. In this 
case, the evolution equations will involve density and magnetization fluctuations allowing for a characterization of 
the phase transitions which lead to the creation of magnetic structures (ferro, antiferro or itinerant) or charge density 
wave instabilities (phase separation or stripe formation) . We believe that a theory able to deal with long wavelength 
fluctuations, like QHRT, may represent a useful tool for tackhng this class of problems. 
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QHRT 


Monte Carlo 


Scri^ 


5"= 1/2 


0.90 


0.94E| 


0.933 


S= 1 


2.66 






5 = 3/2 


5.16 






S = 2 


8.35 






5 = 5/2 


12.25 






S = oo 


1.419 


1.4430 


1.4450 



TABLE I. Critical temperatures for different values of the spin obtained by QHRT and other methods in the three dimensional 
Heisenberg model. 
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FIG. 1. RG trajectories for the two dimensional Heisenberg model computed via numerical integration of the QHRT equation. 
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FIG. 2. Correlation length in two dimensions for S = 1/2. 
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